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Abstract: 

^ ■ The 0{a) electroweak radiative corrections to 77 — > WW — > 4/ within the electroweak 

Standard Model are calculated in double-pole approximation (DPA). Virtual corrections 
are treated in DPA, leading to a classification into factorizable and non-factorizable contri- 
butions, and real-photonic corrections are based on complete lowest-order matrix elements 
for 77 — > 4/+7. Soft and collinear singularities appearing in the virtual and real cor- 
rections are combined alternatively in two different ways, namely by using the dipole 
subtraction method or by applying phase-space slicing. The radiative corrections are 
implemented in a Monte Carlo generator called C0FFER77"'", which optionally includes 
anomalous triple and quartic gauge-boson couplings in addition and performs a convolu- 
tion over realistic spectra of the photon beams. A detailed survey of numerical results 
comprises 0(a) corrections to integrated cross sections as well as to angular, energy, and 
invariant-mass distributions. Particular attention is paid to the issue of collinear-safety 
in the observables. 
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1 Introduction 



As an option at a future e + e~ linear collider, a photon (or 77) collider [1] found 
considerable interest in recent years. It could provide us with information about new 
physics phenomena, such as properties of Higgs bosons or of new particles, which is 
in many respects complementary in the e + e~ and 77 modes (see, e.g., Refs. [1,2] and 
references therein). Moreover, a 77 collider is a true W-boson-pair factory, owing to the 
extremely high W-pair cross section, which tends to a constant of about 80 pb in the 
high-energy limit (in the absence of phase-space cuts), opening the possibility of precision 
studies in the sector of electroweak gauge bosons. Either way, whether one is interested 
in W-boson precision physics or in the search for new phenomena, precise predictions for 
W-pair production are indispensable for signal and background studies. 

In our previous work [3] we have made the first step towards a precision calculation 
for the processes 77 — > WW — > 4/ (+7) by constructing a Monte Carlo event generator 
for lowest-order predictions based on complete matrix elements for the processes 77 — > 4/ 
and 77 — > 4/7. The possibility to convolute the cross sections with realistic photon beam 
spectra is offered upon using the parametrization of CompAZ [4]. The Standard Model 
(SM) predictions were successfully compared to results obtained with the multi-purpose 
packages Whizard [5] and Madgraph [6]. Moreover, we included an effective 77H 
coupling, which is induced by loop diagrams, as well as anomalous triple and quartic 
gauge-boson couplings. The former is needed for studying Higgs production in the s- 
channel. An analysis of anomalous gauge-boson couplings in 77 — > WW provides direct 
information on the 7WW and 77WW interactions without interference from the Z-boson 
sector. Both the Higgs resonance in the s-channel and the more direct access to the 7WW 
and 77WW interactions are complementary to the situation in e + e~ annihilation. 

In this paper we extend our lowest-order calculation [3] for 77 — > 4/ by including 
the electroweak radiative corrections of 0(a) to the W-pair channels 77 — > WW — > 4/ 
in the so-called "double-pole approximation" (DPA). The DPA extracts those contribu- 
tions of the 0(a) corrections that are enhanced by two resonant W-boson propagators, 
i.e. it represents the leading term in an expansion of the cross section about the two 
W-propagator poles. Note that tree-level diagrams for 77 — > 4/ with at most one res- 
onant W boson are suppressed w.r.t. the doubly-resonant 77 — > WW signal by a factor 
of Ofiw/Mw) ~ 0(a). Consequently, predictions based on full lowest-order matrix ele- 
ments for 77 — > 4f and 0(a) corrections for 77 — > WW — > 4/ in DPA should be precise 
up to terms of 0(a/n x Tw/Mw), since corrections typically involve the factor a/n. In- 
cluding a quite conservative numerical safety factor, the relative uncertainty should thus 
be < 0.5% for such predictions, as long as neglected effects are not additionally enhanced. 
The naive error estimate can, in particular, be spoiled by the occurrence of large scale 
ratios, which exist, e.g., near production thresholds or at very high energies. The es- 
timate has recently been confirmed for e + e~ — > WW — > 4/ with centre-of-mass (CM) 
energies 170 GeV < y/s < 300 GeV by comparing a full 0(a) calculation [7,8] with the 
corresponding DPA predictions provided by Racoon WW [9,10]. 

In detail, we apply the DPA only to the virtual corrections to 77 — > WW — > 4/, 
while we base the real-photonic corrections on complete lowest-order matrix elements for 
77 — > 4/7. Apart from the treatment of IR (soft and collinear) singularities, we can use 
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the calculation of the bremsstrahlung processes 77 — > 4/7 for massless fermions described 
in Ref. [3]. The concept of the DPA was already described in Ref. [11] for the corrections 
to e + e~ — > WW — > 4/ and later successfully applied to these processes in different versions 
[9,10,12-14]. We follow the strategy of RacoonWW [9,10] and adapt it to 77 collisions 
where necessary. The virtual corrections in DPA can be naturally split into factorizable 
and non-factorizable contributions. The former comprise the corrections to on-shell W- 
pair production [15-17] 1 and the decay [19] of on-shell W bosons. The latter account for 
soft-photon exchange between the production and decay subprocesses; the known results 
for the non-factorizable corrections [20,21] for e + e _ — > WW — > 4/ can be taken over to 
77 collisions with minor modifications. Although the basic building blocks for the virtual 
corrections exist in the literature, the combination into a complete set of 0(a) corrections 
in DPA has not been done yet for 77 — > WW — > 4/. 

The combination of virtual and real-photonic corrections is non-trivial for two rea- 
sons. First, the finite-fermion-mass effects have to be restored in the phase-space regions 
of collinear photon radiation off charged fermions, and the IR regularization for soft- 
photon emission has to be implemented. To this end, we employ the dipole subtraction 
formalism for photon radiation [22,23] as well as the more conventional phase-space slic- 
ing approach. The second subtlety concerns the fact that we apply the DPA only to the 
virtual corrections, but not to the real-photonic parts. Therefore, the cancellation of soft 
and collinear singularities has to be done carefully, in order to avoid mismatch. 

The paper is organized as follows: After a brief outline of our strategy in the next 
section, in Section 3 we describe the actual calculation of the virtual corrections in DPA. 
Apart from the general concept, we give some details on an efficient way for a numerically 
stable evaluation, on renormalization issues, on the treatment of the s-channel Higgs res- 
onance, and on an improved Born approximation used in the threshold region of W-pair 
production. Section 4 deals with the combination of virtual and real-photonic corrections; 
all relevant details for the application of dipole subtraction and phase-space slicing to the 
considered processes can be found there. Moreover, the differences in the evaluation of 
collinear-safe and non-collinear-safe observables are described. Our discussion of numer- 
ical results is presented in Section 5; besides integrated cross sections, we also discuss 
angular, energy, and invariant-mass distributions. A summary is given in Section 6, and 
the appendices provide further details on the evaluation of coefficient functions for the 
factorizable virtual corrections as well as on the generalization of the dipole subtraction 
method for non-collinear-safe observables. 

2 Strategy of the calculation 

We consider the process 

7 (fc 1 , Ai) + 7 (A; 2 , A 2 ) - W + (A; + , A + ) + W~(A;_, A_) 

-»• h(Pi, o"i ) + h(p2, 02) + h(ps, 0-3) + h(P4, 0-4), (2.1) 

where fcj and denote the momenta and Aj and <Xj the helicities of the corresponding 
particles. 

1 Radiative corrections to on-shell W-pair production, 77 — > WW, were also considered in Ref. [18]. 
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The lowest-order cross section dug^ , based on the complete matrix elements 
•^Bom 4 ^ w hh massless fermions, has been discussed in detail in Ref. [3]. Suppressing 
the averaging over the photon polarizations and the spin and colour summation for the 
final state in the notation, it reads 

/ d^ 4/ = yj ^ f \MH2 4f \ 2 , (2.2) 

with 

s = (h + k 2 )\ s ij = ( Pi + Pj ) 2 , i,j = 1,2,3,4. (2.3) 

The variables are introduced for later use. 

In the following we focus on the radiative corrections of 0(a) which consist of virtual 
corrections dcrj^ 4 ^ to the process (2.1) and real-photonic corrections d<7 77 ~^ 7 , originat- 
ing from the process 

7 (£; 1 , Ai) + 7 (fc 2 , A 2 ) -> W+(A; + , A+) + W~(£;_, A_) ( + 7) 

-»• /l(Pl, + /2(P2, 0-2) + /s(P3, 0-3) + /4(P4, ^4) + l(k, A 7 ). (2.4) 

Combining the different contributions we obtain the C(o;)-corrected prediction for the 
cross section, 

jda = J daH2 4f + / daV^ f + j . (2.5) 

The real-photonic corrections d(T 77 ~ >4 ^ 7 are based on the full lowest-order matrix elements 
-^Born 4 ^ 7 °f the process 77 — > 4/7, which were calculated in Ref. [3] for massless fermions. 
In the limit of vanishing photon momentum k (soft limit) or when the photon becomes 
collinear to an external charged fermion (collinear limit), the cross section diverges. Con- 
sidering the process 77 — > 4/7 with a visible photon (which is neither soft nor collinear), 
these singularities are removed by imposing appropriate phase-space cuts which are jus- 
tified by the finite experimental resolution. For predictions of the 77 — > 4/(7) processes, 
i.e. with or without photon radiation, the singular phase-space regions of soft or collinear 
emission have to be integrated over. In this case the real corrections are combined with 
the virtual corrections which contain exactly the same singularities with opposite sign. 
The regularization of the singularities in the real corrections by small photon and fermion 
masses, A and m,f, as well as the matching with the singularities in the virtual corrections, 
is described in detail in Section 4. The starting point is a separation into a finite and a 
singular part, 

da 77->4/ 7 = + d ^4/ 7; (26) 

where the soft and collinear singularities appear in da^njT 4 ^ 7 as In A and lnmj terms, 
respectively. 

The virtual corrections to the process (2.1) are calculated in the DPA, which is ex- 
plained in Section 3. Since the real corrections are based on complete 77 — > 4/7 matrix 
elements (i.e. they are not calculated in DPA), the cancellation of soft and collinear sin- 
gularities in Eq. (2.5) requires particular care. To this end, we apply the DPA only to the 
finite part of the virtual corrections, 

uo virt ^ ut virt,fimte,DPA ^ ao virt,sing- l z - ' J 



3 



Technically this is achieved by subtracting the singular part in DPA from the DPA virtual 
corrections and adding the exact singular part da"^~^f g . Of course, this procedure involves 
some freedom, because finite terms can be shifted between du^~^ teDPA and da^rt^ng- 
This arbitrariness is, however, of the order of the uncertainty 0{aY^ / {nM^)) of our 
calculation. In the e + e~ case this has been checked numerically in Ref. [10]. 
Inserting these rearrangements into Eq. (2.5) we obtain 

J = J dfTBorn ^ + J d^virt.firlrte.DPA + J ^^virt+real.sing + J ^Hite ^ > (2-8) 

where / da^ 4 r ii, S in g = / daV^l g + 1 da rYai2S does not contain any dependence on the 
photon mass anymore. Collinear singularities, appearing as mm/ terms, also cancel if 
the observable is sufficiently inclusive. Such collinear- safe observables result if photons 
within cones collinear to any outgoing charged fermion are treated inclusively, i.e. if they 
are not separated from the nearly collinear fermion by any phase-space or event selection 
cuts. For non-collinear-safe observables logarithms of the fermion masses remain in the 
final result. This case demands a special treatment of the singular terms. We elaborate 
more on this issue in Section 4.2. 

The radiative corrections are implemented in a Monte Carlo generator called 
C0FFER77, which is based on the lowest-order calculation described in Ref. [3]. We 
emphasize that we have actually constructed two independent Monte Carlo programs, 
each of which employs independent routines for the matrix elements (with and without 
corrections), for the subtraction procedure, and for the phase-space integration. The nu- 
merical results obtained with the two programs are in mutual agreement within statistical 
uncertainties. 

3 Virtual corrections 

3.1 Concept of the double-pole approximation 

In the DPA the matrix element for 77 — > 4/ is expanded around the poles of the two 
resonant W propagators. The leading term of this expansion receives contributions from 
so-called factorizable and non-factorizable corrections. For the details of this classification, 
especially how a gauge-invariant decomposition is obtained, we refer to Refs. [10,11,20,21]. 

The generic Feynman diagram for the factorizable corrections is shown in Figure 1. 
It factorizes into the on-shell W-pair production, the off-shell W-boson propagators, and 
the subsequent on-shell W decays. The corrections can be attributed to either of these 
subprocesses. When integrating over the full 4/ phase space, the W bosons usually are not 
on shell. However, a gauge-independent evaluation of the matrix elements for production 
and decay requires on-shell momenta for the W bosons. Therefore, we have to perform an 
on-shell projection, i.e. the momenta of the fermions are deformed in such a way that the 
W bosons become on shell. The deformation involves a certain freedom and introduces 
an error of (9(ar w /(7rM w )). We define the on-shell projection by fixing the directions 
of the W + boson and of the fermions fi and / 3 . The explicit formulas can be found in 
Appendix A of Ref. [10]. For later use, we label the new momenta k± and pi and define 
the kinematic invariants 

t=(k 1 -k + ) 2 = (k 1 -p 1 -p 2 ) 2 , u = 2M^-s-l (3.1) 
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On-shell production On-shell decays 



Figure 1: Generic Feynman diagram of the virtual factorizable corrections to 77 — > 
WW — > 4f. The shaded blobs stand for loop corrections to the production and decay 
processes. 

Apart from the factorizable corrections there are additional doubly-resonant contribu- 
tions. In the corresponding diagrams subprocesses are linked by a photon. These diagrams 
become doubly resonant in the limit of vanishing photon momentum, as can be seen from 
the soft-photon approximation in which the correction is proportional to the lowest-order 
cross section. The relative correction factor for these so-called non-factorizable correc- 
tions is, thus, not dependent of the actual production mechanism of the W pairs, but 
only on the electric charges and kinematics of the external particles of the process. The 
non-factorizable corrections were calculated in Refs. [20,21] for e + e~ — > WW — > 4f. We 
can transfer the results for the e + e~ case by simply omitting all contributions in which 
the exchanged photon is linked to an from the initial state. The different types of 
relevant diagrams are depicted in Figure 2. The first two diagrams, labelled (a) and (b), 
are manifestly non-factorizable, i.e. the photon links different subprocesses so that the 
propagators in the diagrams cannot be factorized anymore. The diagrams (c), (d), and 
(e) contain both factorizable and non-factorizable contributions. Their factorizable parts 
are defined as the residues for on-shell W bosons times the off-shell W-boson propaga- 
tors; note that this procedure introduces artificial soft IR divergences connected with the 
on-shellness of the W bosons in the loops. The non-factorizable parts of the diagrams are 
obtained from the difference of the doubly-resonant contribution of the full diagrams and 
their factorizable parts; the artificially introduced IR divergences of the factorizable parts 
are, thus, compensated by corresponding terms in the non-factorizable parts. 

Following this strategy, the virtual corrections in DPA can be written as 

d^iIt?DPA^ 4/ = d$ 4/^2Re{(5A4 v irt,f ac tA4Born,DPA} 

+ "virt,nfact \M Born.DPA r+i^mggsi 2 ), (3.2) 

where A^Bom,DPA denotes the tree-level matrix element in DPA and <5 v irt,nfact contains 
the non-factorizable corrections. The factorizable corrections 5M. vlvt ^ act also contain a 
contribution of the s-channel Higgs resonance, ^A^Higgs- In order to describe this resonance 
properly, it is not sufficient to include the interference of SAi Hi ggs with the Born matrix 
element, but the square of this matrix-element contribution has to be taken into account in 
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addition. To this end, dMmggs has to be defined in a gauge-invariant way. Our treatment 
of ^.Mffiggs is described in Section 3.2.4 in detail. 

3.2 Factorizable corrections 

3.2.1 Calculation of the one-loop amplitudes 

The factorizable corrections comprise the corrections to the on-shell production of the 
W bosons and their on-shell decay and can be expressed as 

JAW* = E j^- {8M^-MZ7 hh MZ- hh 

+ MlZ +W ~SM w+ ^MZrr hf4 

+ MlZ +W ~M^ hf2 SM^-^) , (3.3) 

where we introduced the abbreviations 



K± = k?±-M%, + iM w r w , (3.4) 

and SAi denote one-loop matrix elements. Note that all matrix elements on the r.h.s. 
of Eq. (3.3) depend on the on-shell projected momenta, but the momenta in K± remain 
unchanged. The results for the different one-loop corrections are already known in the 
literature [15-17,19]. Combining them in Eq. (3.3) is, however, non-trivial since the 
polarizations of the W bosons have to be defined consistently in a common reference 
frame. 

The one-loop corrections 5.M W_> -^ to the W decays are rather simple. In the massless 
limit they are proportional to the respective Born matrix elements Ai^^f 1 * 3 , 

6M^*H><vr,Pi,Pi) = ^ fJj MZ7J Jj ^w,Pi,Pj), (3.5) 

where <5 W ~*-^ is a constant correction factor that neither depends on the kinematics nor 
on the helicity Aw of the decaying W boson. 

The one-loop correction <5.A/1 77W+W to the W-pair production process contains the 
complicated part; we have derived it in two independent ways. One calculation is based on 
the results of Ref. [16] which were obtained in a gauge with a non-linear gauge fixing term. 
The other is a new calculation based on the program FeynArts [24] for the generation of 
the amplitudes and on in-house Mathematica routines for their algebraic reduction. This 
second calculation has been carried out in the 't Hooft-Feynman gauge and repeated 
in the background-field gauge [25] to get an additional consistency check. The results 
obtained from the different calculations are in mutual numerical agreement. 

In the following we describe an efficient way for calculating the contribution of 
SAi llW+w ~ to £.M v i rti fact of Eq. (3.3), taking into account all spin correlations. As de- 
scribed in Ref. [16], the matrix element <5.A/f 77W+w for on-shell W-pair production is 
decomposed into a sum of products of form factors Fj, which only depend on the kine- 
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matic variables s and i, and a set of standard matrix elements (SME) jny w+w ) which 
contain the polarizations and momenta of the external photons and W bosons, 

5M 7JW+W- ^ fc2> £ +> ^ A+; A ) 

36 

= £ ^(s, t)-Mj 7W+W_ (fci, fc 2 , A 1; A 2 ; k + , fc_, A + , A_). (3.6) 
i=i 

The SME _A^7 7W+W are obtained from the 83 basic matrix elements given in Section 2 
of Ref. [16] which are reduced to 36 matrix elements as described there 2 . The decay 
matrix elements -M^,„ fj , which multiply 5.M 77W+W in Eq. (3.3), can be included by 
replacing the W polarization vectors e*± in the definitions of the SME Aij by the "effective 
polarization vectors" 

4" = ^^u(Pih^-v(P2), tJ = '-LzfaWu.vfa), (3.7) 

where u(pi) and v(pi) are the Dirac spinors of the fermions and antifermions and uo_ = 
|(l—75) is the left-handed chirality projector. The effective W-polarization vectors e*± are 
formal shorthands for the W propagators and the tree-level decay matrix elements, which 
involve the usual SU(2) gauge coupling e/s w . Upon substituting e* ± — > s*± in the SME for 
on-shell W-pair production, we obtain a new set of SME Mj that correctly transfer the 
W polarization to the decay, 

Mjih, k 2 , Ai, A 2 ; k\ , k 2 _\ {ft}) = .MJ 7 ^™"^!, k 2 , X u A 2 ; fc+, fc_, A+, A_) 

x Mtrn hf ~ 2 (A+, pi, ft) ACrT^ 4 (A-, ft, p 4 ) ■ 

(3.8) 

The new SME yVfj can be easily evaluated with spinor methods, as e.g. described in 
Ref. [26]. 

In summary the factorizable part of the virtual correction takes the form 

36 

<5-M vir t,f ac t = E F j( s > t)Mj{h, k 2 , Ai, A 2 ; k 2 + , k 2 _; {pi}) 
i=i 

+ (sw+^fj, + s w-^f 3 f^ M Boin , DPA (k u k 2 , A 1? A 2 ; A; 2 , P ; {p,}). (3.9) 

3. 2.2 Details of the numerical evaluation 

The formulas for the coefficient functions Fj are rather lengthy and contain many one- 
loop integrals, which in turn involve many dilogarithmic functions, etc. Thus, to speed 
up the numerical evaluation it is desirable not to evaluate the Fj at each phase-space 
point. Moreover, numerical instabilities occur at the boundary of the phase space where 



2 The on-shell momenta k± and the helicities A± of the W bosons are denoted £3,4 and A3, 4 in Ref. [16]. 
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the scattering angle 9 between the W bosons and the beam axis tends to or n. This is 
due to the inverse Gram determinants appearing in the Passarino-Veltman reduction [27] 
of the tensor integrals. The problems of speed and stability can be solved by expanding 
the functions Fj(s,t) in terms of a generalized Fourier series in the variable t for fixed 
values of s. The coefficients of this expansion are calculated before the Monte Carlo 
integration. An appropriate system of orthogonal functions in the variable x = cos9, 
which is equivalent to a function of i for fixed s, is provided by the Legendre polynomials 

Pl{x) = \- {x " - 1} 1 ' l = ^- ■ ( 3 - 10 ) 
For this basis functions, the coefficients read 

21 + 1 r+ 1 

cM = —^~ J dcos9(t-M^)(u-M^)F j (s,t)P l (cos9), (3.11) 

where we have introduced the factor (t — M^){u — M^) in order to flatten the t- and 
w-channel poles in the functions Fj. This improves the efficiency of the expansion. The 
integration in Eq. (3.11) is carried out using Gaussian integration. With 40 integration 
points the region of instability is not entered (for energies up to a few TeV), and the 
integration is sufficiently precise. During the Monte Carlo integration the coefficient 
functions are recovered by the generalized Fourier series 

OO 1 

Fj(s,t) =Y-~ 7— 7- c,- , (s)P, (cos 0). (3.12) 

In Ref. [10] the same concept was used to evaluate the factorizable corrections to 
e + e~ — > WW — > 4/; there it was sufficient to use the Legendre polynomials up to / = 20 
for a good accuracy. In the case of 77 — > WW, however, the coefficient functions involve 
inverse Gram determinants l/(tu — M^) oc 1/ sin 2 9 which appear in the Passarino- 
Veltman reduction of the tensor integrals. As each step in this recursive reduction involves 
such an inverse determinant, l/(tu — M^) can appear up to the fourth power. At cos 9 ss 
±1 this factor leads to a behaviour of the Fj(s,t) that is not well approximated by the 
Legendre expansion. Using higher-order Legendre polynomials is not a solution since this 
increases the calculation time and also requires more integration points for the Gaussian 
integration. The more points are used in the Gaussian integration, the closer some of these 
points approach the integration boundary where the numerical stability of the coefficient 
function breaks down. Therefore, we follow a different strategy based on the fact that the 
helicity amplitudes for the on-shell process 77 — * WW are smooth functions of cos 9, apart 
from the t- and w-channel poles. Thus, within the full amplitude the factors l/(tu — M^) 
have to cancel between contributions of different coefficient functions. To make use of this 
fact we change the basis of SME by a linear transformation in such a way that the new 
coefficient functions correspond to helicity amplitudes of the on-shell process 77 — > WW. 
Some details of this transformation can be found in App. A. After this transformation 
the uncertainty of the approximated matrix elements in Eq. (3.6) is well below 10~ 4 with 
respect to the Born matrix elements for all values of cos 9. 

In contrast to the e + e~ case, the CM energy ^Js of the photons is not fixed. Thus, we 
have to perform the Legendre expansions for different values of s. During the Monte Carlo 
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integration we derive an approximate value of the coefficients Cjj{s) by interpolation. Since 
the Fj(s,t) depend on s very smoothly, it is sufficient to calculate the Cjj(s) at intervals 
of As < 1 GeV. In these intervals we then interpolate with a polynomial of third order. 
We have checked that, up to 1 TeV, this yields a sufficient accuracy (i.e. better than the 
accuracy of the Legendre expansion). 



3.2.3 Renormalization and imaginary parts of virtual corrections 

For on-shell W-pair production, which was considered in Ref. [16], imaginary parts of 
counterterms, if included, do not influence the correction to the matrix element square. 
The reason is that for the 2 — > 2 scattering process 77 — > WW all SME, and thus also the 
Born matrix element, can be taken real by appropriate phase choices. Thus, the operation 
of taking the real part in the interference term 2 Re{A^ ct A^B orn } of the counterterm 
contribution M. ct to the one-loop amplitude with the Born amplitude effectively acts on 
the renormalization constants themselves. The same argument shows that also imaginary 
parts of loop integrals drop out. These arguments are no longer true if the decay of the 
W bosons is taken into account, because the SME and the Born matrix element A^Bom.DPA 
become necessarily complex. Thus, imaginary parts of renormalization constants and of 
loop integrals in general matter. Considering the W-decay amplitudes in the DPA in 
more detail, as e.g. done in Ref. [13] for the e + e~ case, one can see that imaginary parts 
average to zero after the azimuthal decay angles of the W-decay products are integrated 
over. 

We have calculated the virtual corrections taking into account the imaginary parts of 
all loop integrals. As already mentioned, we carried out the whole loop calculation in 
different gauges: in the 't Hooft-Feynman gauge [28], in a non-linear gauge [16], and in 
the background-field gauge [25] . We find agreement between the results obtained in these 
different gauges, but only if we also take into account the imaginary parts of the loops 
that contribute to renormalization constants. In order to explain this fact, we consider 
the counterterm contributions to the one-loop matrix element in more detail. 

Following Ref. [16], we write the Born matrix element in DPA as 

A^Bom,DPA = 8na \ — - M ,t + — -M , u - (eie 2 )(e^e*) \ , (3.13) 

where A4o,t and A4o,u are abbreviations for specific combinations of momenta and po- 
larization vectors defined as in Eq. (22) of Ref. [16] for on-shell W-pair production. In 
the 't Hooft-Feynman gauge, the counterterm contribution to the production part of the 
factorizable correction reads 



^^ctlprod — A^Born.DPA [25Z e + 5Z W + 5Z AA -5Z ZA \ 

\ s w J 



3 w 

o ( sSM w 1 j sSM^ N 

V (t - (u - M$f) J V SwC w J 

+ 4 ™^7 {(i-M^ + (u-M^y) 5t > (3 - 14) 
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where we adopt the conventions of Ref. [28] for the renormalization constants SZ e , SZ W , 
etc. The explicit calculation of the constants in terms of self-energies is also described 
there. The counterterm contribution in the background-field gauge [25] can be obtained 
from 5-A/f*^p rod by simply omitting the SZza terms, because 5Zza vanishes owing to the 
background-field gauge invariance. In the non-linear gauge the counterterm contribution 
reads 



as described in Ref. [16], which is different from its counterpart in 't Hooft-Feynman 
gauge. Note also that the explicit expressions of the renormalization constants in the 
different gauges are in general different. 

Imaginary parts of loop and counterterm contributions that are proportional to the 
Born matrix element, SM. = cA^Born, cannot influence matrix element squares, because 
2 Re{SMM% OTn } = 2 Re{c}|.MBorn| 2 - Thus, the W-mass renormalization constant SM^ 
is the only renormalization constant whose imaginary part plays a role, since the tadpole 
counterterm St is a real quantity. From Eqs. (3.14) and (3.15), we see that SM^, which is 
equal in all three considered gauges, enters the counterterm contributions in the 't Hooft- 
Feynman gauge and in the non-linear gauge in different ways. In fact, we have checked 
numerically that the virtual corrections in these two gauges are different (though finite) 
if the usual on-shell prescription SM^ = Re{Y^ (M^)} (see e.g. Ref. [28]) is applied, 
where T^(k 2 ) is the transverse part of the W-boson self-energy with momentum transfer 
k. If we, on the other hand, use the definition SM^ = EJjf(M{y), i.e. without taking the 
real part of the self-energy, we find agreement for the results from the different gauges. 
This clearly shows that the imaginary part of a one-loop amplitude is in general gauge 
dependent if imaginary parts in renormalization constants are not taken into account. 
The reason for this fact, in other words, is that the decomposition of a renormalized 
transition matrix element into genuine loop parts and counterterm contributions depends 
on the gauge fixing. 3 

In our Monte Carlo generator we have taken into account the imaginary parts of the 
virtual corrections (including the ones from counterterms) ; more precisely they can be 
switched on and off optionally. As explained above, they could only affect observables 
that are sensitive to the azimuthal decay angles of the fermions. In our numerical results, 
we could, however, find no significant effects. 

3.2.4 Higgs resonance 

The loop-induced Higgs resonance, 77 — > H — > WW — > 4/, belongs to the class of 
factorizable contributions. Nevertheless, its treatment, especially the question of gauge 
invariance when including the Higgs decay width, deserves some care. In Ref. [16] the 

3 A consistent renormalization prescription with complex renormalization constants naturally leads to 
complex masses for unstable particles. Such a renormalization scheme was proposed recently in Ref. [8] 
in the context of a full 0(a) calculation for e + e _ — > 4/. 



SM 



ct,prod 



NL 




(3.15) 
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^Hi gg s = — „ 2 , ;il/f ^ • (3.19) 



diagrams with an s-channel Higgs resonance were decomposed into a gauge-invariant 
resonant part and a gauge- dependent non-resonant part. If we write the contribution of 
the Higgs-exchange diagrams as 

5M~^ = I^( £l e 2 )(iia, (3.16) 

with F H (s) given in Section 4.3 of Ref. [16], and e 1 and e 2 being the polarization vectors 
of the photons, then the Higgs decay width can be introduced by replacing 

^ y s — M\ + iMnTn ^Mg ) ^e 2 ){e + e_). (3.17) 

As the residue F H (M^) is gauge independent, we have introduced the Higgs decay width 
Th in a gauge-invariant way. Recall that the choice of the polarization vectors of the 
photons is such that they obey 

5.^ = 0, i,j = 1,2. (3.18) 

Close to the resonance, the contribution of the Higgs-exchange diagrams is strongly 
enhanced. This is why we also take into account the square of the resonant part in 
Eq. (3.2), 

F*(Mg)( £l£2 )(e;e*_ 
s - Ml + iM H r H 

In this approach only the leading contribution to the Higgs resonance is taken into 
account. However, the gauge- invariant separation of SA4 Hi ggs from the remaining one- loop 
amplitude easily allows for specific improvements in predictions for the Higgs-production 
signal in the future. To this end, a pole expansion about the Higgs resonance would 
be an adequate first step. Conceptually this expansion again leads to factorizable and 
non-factorizable contributions, but the corresponding ingredients are not all available yet 
and their calculation is beyond the scope of this work. It should be mentioned that both 
the 0(a) electroweak and 0{a s ) QCD virtual factorizable corrections to (on-shell) Higgs 
production 77 — > H can be deduced from the corresponding two-loop calculations [29] 
(see also references therein) for the decay H — > 77. 

3.3 Non-factorizable corrections 

As explained in Section 3.1, we make use of the result for the non-factorizable correc- 
tions to e + e~ — > WW — > Af. According to Refs. [10,21] we write the correction factor 
to the lowest-order cross section as a sum over contributions that are associated with 
different pairs of fermions, 

<W&ct= E Y.(- l ) a+b+1 QaQb-^{^\k + ,p a] k-, Pb )}. (3.20) 

a=l,2 6=3,4 71 

The function A virt receives contributions from the different types of diagrams in Figure 2, 

A vi rt = A virt + A virt + A virt + A virt + ^ (3^) 
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(a) type (mf' 




(b) type (ff) 




(c) type (mm') 




(d) type (mf) 




(e) type (mm) 




Figure 2: A representative set of diagrams contributing to the virtual non-factorizable 
corrections. The shaded blobs stand for all tree-level structures contributing to 77 — > 
WW. 

for which the results were given in terms of scalar integrals in Ref. [10]. The final result 
for a = 2, b = 3 (all other contributions can be derived by appropriate substitutions) is 



Avirt 1 a virt , A virt 
mf + A ff + A mf 



A 



virt 
mm' 



D {-P4, k + + p 3 , p 2 + p 3 , 0, M, M, 0) 



K + K_s 23 det(Y 
det(y) 



det(y) 



det(Y) 



S23 . \ 



(2M^ - s) C (k + , -k., 0, M, M) - C (k + , -k_, A, M w , M w ) 



k 2 ± =M w 



2 [ ' 



21n^ + 21n^ + 4, 



—K. 



(3.22) 
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where the sign "~" indicates that the limit k\ — > and T w — > is carried out whenever 
this does not lead to a singularity. The matrices Y , Y 2 , Y 3 , and Y arise from the reduction 
of 5-point functions and can be found in Section 3.1 of Ref. [21]. The functions F 2 and 
F 3 are defined in Section 4.2, and the Co and D functions in Appendix C.l of the same 
reference. The contribution A^, contains the difference of the full off-shell and on-shell 
Coulomb singularity, as described there in detail. 

The full correction factor 5 v i r t, n fact does not contain fermion-mass singularities [10], but 
involves IR-singular terms In A, as explicitly visible in Eq. (3.22). The latter originate from 
the subtraction of the virtual factorizable correction, which involves the one-loop matrix 
elements for 77 — > WW and W — > //' with on-shell W bosons, from the doubly-resonant 
part of the matrix element for the full 77 — > 4/ process. Specifically, the In A terms stem 
from diagrams with photon exchange between an on-shell W boson and another on-shell 
particle. As already explained in Section 3.1, these singularities cancel in the sum of 
factorizable and non-factorizable contributions, since they are artificially introduced in 
the corresponding decomposition of the virtual correction. 

3.4 Leading universal corrections and input-parameter scheme 

We parametrize the cross section in such a way that the universal corrections arising 
from the running of the electromagnetic coupling a and from the p-parameter are absorbed 
in the lowest order. To this end, we take all particle masses as input, from which the 
weak mixing angle is derived via the on-shell condition 

^ = 1 -4 = 1 -§- (3 - 23 > 

The electromagnetic coupling a is chosen in order to absorb some universal corrections. 

As pointed out in Ref. [16], the relevant coupling for the 77 — > WW production 
process is the fine-structure constant a(0), because the external on-shell photons do not 
induce any running in their coupling to the W bosons. This means that the on-shell 
renormalization is carried out precisely as described in Refs. [25,28] for this contribution. 

For the decay of the W bosons, it is, however, appropriate to derive a from the Fermi 
constant leading to 

a G = ^ G " M w4 (3 24) 

M 7T 

This modification of the coupling induces an additional finite contribution to the charge 
renormalization constant, 

5Z C \ G ^ = 5Z c \ a{0) - -Ar, (3.25) 

where 5Z e \ a ^ is the charge renormalization constant of the on-shell renormalization 
schemes [25,28] with a(0) as renormalized coupling. The quantity Ar contains the ra- 
diative corrections to muon decay; explicit expressions for Ar can, e.g., be found in 
Refs. [28,30]. 

In summary, our lowest-order cross section scales like a(0) 2 a;^. For the relative 0(a) 
corrections we use a(0), which is the correct effective coupling for real photon emission, so 
that the corrected cross section scales like a(0) 3 a^. For the loop-induced Higgs resonance 
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we exceptionally take the scaling factor a(0) 2 o;^, which accounts for the two "photonic" 
and the three "weak" couplings in the corresponding diagrams. We perform this rescaling, 
of course, only in the gauge-invariant resonant part <5.Mffi g gs of the one- loop amplitude, 
as defined in Eq. (3.19). 

3.5 Improved Born approximation 

The motivation for calculating the virtual corrections in DPA lies in the domi- 
nance of doubly-resonant diagrams. At threshold, however, singly-resonant and non- 
resonant diagrams become equally important, thus, rendering the naive error estimate of 
0(aiTyf / (My/ir)) unreliable. As a consequence, we decided to use the DPA only for a CM 
energy v /s^ 7 " > 170 GeV when integrating over the photon spectrum. For ^/s 77 < 170 GeV 
we make use of an improved Born approximation (IBA), i.e. we approximate the 0(a) 
corrections by universal contributions without any expansion about the W resonances. 
Assuming that the IBA accounts for all 0(a) corrections with pronounced enhancement 
factors, the relative uncertainty of the IBA is about ~ ±2%. For the corresponding 
e + e~ reaction this expectation was confirmed by the full 0(a) calculation [7,8] for 4/ 
production. 

In detail, we start from the Born cross section based on the full set of 77 — > 4/ 
diagrams, which is parametrized as described in the previous section and include the 
Higgs resonance with SM couplings. Such lowest-order predictions, which we denote 
"Born+Higgs" below, have already been presented in Ref. [3]. In addition, we now dress 
the resulting cross section with the off-shell Coulomb singularity, 

/ <KbT 4/ = y s I d$ 4/ (1 + 5 coul )\Mll2% igg f- (3-26) 
The correction factor <5 cou i for the Coulomb singularity was calculated in Ref. [31] to 

with the abbreviations 

P = -^s 2 + (k\) 2 + (k 2 _ ) 2 - 2sk 2 + - 2sk 2 _ - 2k 2 + k 2 _, 



1 _4(M^-iM w r w); A= |*e-*fl (3 . 28) 
V s s 

4 Treatment of soft and collinear photon emission 

We calculate the real photonic corrections from the full lowest-order matrix element 
of the process (2.4) without any expansion about the W-boson resonances. They are 
calculated from the integral 

J do -77-4/ 7 = J_ J 1^77-4/^2 Q ($ 4/7 ), (4.1) 

where we have made the implementation of phase-space cuts explicit by including the step 
function 0($ 4 / 7 ), which is equal to 1 if an event passes the cuts and otherwise. Since 
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we evaluate the real matrix element yVf 77 ^ 4 / 7 w jth massless particles, the phase-space 
integral diverges in the soft and collinear regions, where the emitted photon is either soft 
or collinear to an outgoing external charged fermion. In these regions we reintroduce a 
formally infinitesimal photon mass A and small fermion masses mj as regulators. 

To this end, we apply two different methods: the dipole subtraction and the (two- 
cutoff) phase-space slicing methods. In the case of collinear-safe observables we closely 
follow the approach of Ref. [10] and only give a brief description in Section 4.1 since the 
procedure is very similar to the e + e~ case. In Section 4.2 we describe how the two methods 
are extended to non-collinear-safe observables. 



4.1 Collinear-safe observables 

4-1.1 Phase-space slicing 

In the phase-space slicing approach the phase space is divided into regions where 
the integrand is finite and can, thus, be integrated numerically, and regions where the 
integrand becomes singular. In the singular regions the integration over the photon phase 
space is carried out analytically in the approximation that the photon is soft and/or 
collinear to a charged fermion. 

The singular regions consist of two parts one of which contains a soft photon (k° < AE) 
and the other a photon that is collinear but not soft (A; > AE and 9 1 f < AO, where 1 j 
is the angle between the photon and a charged fermion). Thus, the real corrections are 
decomposed according to 



J d(T 77-4/ 7 = J d(Jsoft + J d(7coU + J 



, 77^4/7 
uo finitc ) 



(4.2) 



where the cutoff parameters AE and AO are defined in the CM system of the incoming 
photons. Both in the soft and collinear regions the squared matrix element |vVi 77 ^ 4 ^ 7 | 2 
factorizes into the squared lowest-order matrix element \M-'b ( ^ l ^\ 2 and a universal factor 
containing the singularity. The five-particle phase space also factorizes into a four-particle 
phase space and a photon part, so that dcr so ft and dcr co u can be integrated over the photon 
momentum. Taking over the results from Ref. [10] yields 



cKoft = daH2 4f 6($ 4/ ) *±± (-l^QiQj | 2 In (^) 

11 i=l j=i+l { \ A / 



In 



4p°p°' 



In 



ir 3 



miirij 



+ ln 2 



rrii 



&ij 



(4.3) 



and 



d^coii = d^ 4/ e($ 4/ ) — YM 



3 , /A£\ 
-+21n — =- 
2 



l-21n 



'AQp® 



rrii 



2tt 5 



(4.4) 



where Qi and rrii denote the relative electric charge and mass of fermion /j, respectively. 
The step function 0($4/) indicates that both do" so f t and d<7 co ii are defined on the four- 
particle phase space of the lowest-order cross section, so that the singular part 



dcr^g 



.4/7 



da- so ft + ddcoll 



(4.5) 
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can be locally combined with the singular part of the virtual corrections, which are defined 
on the same phase space. In the result do"^j~^{ al sing all dependences on the photon and 
fermion masses cancel. 

While du^^eai sing depends on the cutoff parameters AE and AO analytically, the 
finite real corrections J dcrg^ 4 ^ 7 only show this dependence upon the cuts in the numerical 
integration. Nevertheless, the cutoff dependence has to cancel in the full result in the 
limit AE, AO — > 0. This is illustrated on the l.h.s. of Figures 3 and 4 where the relative 
correction factor 5 = a/a Born - 1 of the 4/ part (/ dCT^ fi 4 n { teDPA + / d<7 rt + 4 r { aljSing ) and 

of the 4/7 part / dcTg^ 4 ^ 7 is shown as a function of the cutoff parameters AE and 
A0. The cancellations of the cutoff dependence of the two contributions is shown on a 
smaller scale on the r.h.s. of Figures 3 and 4. While terms of 0(AE/E hcSLm ) and O(A0) 
become visible for large values of the cutoff parameters, for smaller values a plateau is 
reached. The integration error increases with decreasing cutoff values, until for too small 
values the integration error is usually underestimated. As a result, we decided to take 
AE/E hcam = KT 3 and A6 = 1(T 2 as default values. 



4-1.2 Dipole subtraction method 

In a subtraction method an auxiliary function is constructed that contains the same 
singularities as the real corrections. Subtracting this function from the real corrections, 
this difference can be integrated numerically. The next step is to perform the singular 
integration of the auxiliary function over the photon momentum analytically and to readd 
the result to the virtual corrections. In our case where soft and collinear singularities 
originate from final-state radiation only, the soft and collinear singularities completely 
cancel against their counterparts in the virtual corrections for collinear-safe observables. 

In the dipole subtraction method [22,23], which was originally proposed for QCD [32], 
the auxiliary function consists of different contributions labelled by all ordered combi- 
nations of two charged fermions % and j, which are called emitter and spectator. These 
contributions contain the singularities connected with the emitter %. Since there are 
only charged particles in the final state in 77 — > 4/, the situation is simpler than for 
e + e~ — > 4/. Explicitly the auxiliary function, which is subtracted from the spin-summed 
squared bremsstrahlung matrix element, reads 

4 

l-Msubl = X! I-Msub,y| , 



*>j=i 



|AW*4/7)I 2 = -(-V i+j QiQiZ 2 9T\Pi^ k)\MllZ if (^)\ 2 - (4.6) 

Adopting the formulation of Ref. [22] 4 , the soft and collinear divergences are contained in 
the function 



(pik)(l-yij) 



z ij{}- Uij) 



Zij 



(4.7) 



4 Thc formulation of Rcf. [23] differs from that by the regular factor 1/(1 — y^) in Eq. (4.7), so that 
the readded singular contributions of Refs. [22] and [23] differ by non-singular finite parts. 
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Figure 3: Dependence of the corrections on the energy cutoff in the slicing approach for 
the process 77 — > z/ e e + du at ^s 77 = 500 GeV. For comparison the corresponding result 
obtained with the dipole subtraction method is shown as a la band in the plot on the 
r.h.s. 
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Figure 4: Dependence of the corrections on the angular cutoff in the slicing approach for 
the process 77 — > z/ e e + du at ^/s 77 = 500 GeV. For comparison the corresponding result 
obtained with the subtraction method is shown as a la band in the plot on the r.h.s. 
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with 

V* = V A r> % = -^-r- (4-8) 

PiPj + Pik + Pjk pipj +pjk 

The embedding of the 4/ phase space $4/,ij into the 4/7 phase space $ 4 / 7 is defined as 



with all other momenta unchanged, = Pk,k i,j. Subtracting the auxiliary function 
from the real corrections enables us to carry out the numerical integration, 



(4.10) 



which does not contain any soft or collinear divergences by construction of |-M su t,| 2 for 
collinear-safe observables. In this context, it is important to notice the different arguments 
of the step functions 9 which account for phase-space cuts. Since for a generic point in 
4/7 phase space each ij contribution corresponds to a different point in phase space, there 
is in general no correlation between the values of the different step functions. For collinear- 
safe observables, however, we have 0($ 4 j 7 ) = 6($4/,ij) in the soft region (k — > 0) and in 
the region where the photon momentum k is nearly collinear to the emitter momentum 
Pi (pik — > 0). The collinear safety can, e.g., be enforced by photon recombination, as 
discussed in the next section in more detail. 

In order to combine the subtraction function with the virtual correction, it has to be 
integrated over the photon momentum, yielding 

/ daV- 4f " = -|- £ (-lf^QiQ- J d$ 4/ 4 sub) ( % )|A!Srn 4/ ($4/)| 2 e(<I. 4/ ). (4.11) 

The singularities are contained in the function 

G ( t\s tl )=C(s ll ,ml)-^ + l (4.12) 

with 

£(8y , m?) . ln g) ln g) +ln g)-i 1 ,g) + I ln g). (4 ,3) 

We have numerically checked that these soft and collinear divergences are completely 
cancelled by their counterparts in the virtual correction. 

4.2 Non-collinear-safe observables 

In the previous sections the matching of real and virtual corrections was described 
for collinear-safe observables. We speak of collinear-safe observables if a nearly collinear 
system of a charged fermion and a photon is treated inclusively, i.e. if phase-space selection 
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cuts (or histogram bins of distributions) depend only on the sum p { + k of the nearly 
collinear fermion and photon momenta. In this case the energy fraction 



p\ + k° 



(4.14) 



of a charged fermion /j after emitting a photon in a sufficiently small cone around its 
direction of flight is fully integrated over, because it is not constrained by any phase-space 
cut (or histogram bin selection in distributions). Thus, the KLN theorem [33] guarantees 
that all singularities connected with final-state radiation cancel between the virtual and 
real corrections, even though they are defined on different phase spaces. A sufficient 
inclusiveness is, e.g., achieved by the photon recombination described in Section 5.1, which 
treats outgoing charged fermions and photons as one quasi-particle if they are very close 
in angle. 

In the previous section we could, therefore, integrate the subtraction function |.M su b| 2 
and the slicing contribution d<j co n over Zi analytically. In this section we are concerned with 
non-collinear-safe observables, i.e. the fermion-photon system is not treated inclusively 
and fermion-mass singularities can become visible. As the integration over Zi now is 
constrained by phase-space cuts (or histogram bins), we have to modify the methods 
described in the previous section in such a way that the integration over is part of the 
numerical phase-space integration. 

4-2.1 Phase-space slicing 

In the slicing method the procedure is straightforward. The numerical integration over 
z = Zi in the collinear parts reads 



2tt i=1 jo 



1-AE/pO 



dz &(pi = zpi, k = (1 - z)pi, {Pj^i}) 



x \P f f(z) 



2 In ^z 



rrii 



+ (1-*) , 



with the splitting function 



1 + z 2 



(4.15) 



(4.16) 



The Born cross section and the logarithm still depend on the momenta of the 4/ phase 
space $4/ which are labelled pi- In the cut and recombination function 6, however, the 
momentum pi of the fermion % (before photon emission) is distributed to the fermion mo- 
mentum pi and the photon momentum k. For collinear-safe observables, as e.g. achieved 
by photon recombination, the 6 function effectively only depends on the sum Pi + k = pi 
of the collinear momenta, which is independent of z. In this case, the 6 function becomes 
0($4/), and the ^-integration can be easily carried out analytically yielding Eq. (4.4). 



4-2.2 Dipole subtraction method 

In the case of the dipole subtraction method the generalization to non-collinear-safe 
observables is more complicated than in the slicing approach, since the integration over 
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the photon momentum is more involved. Here, we collect the formulas relevant for our 
calculation. Details on their derivation are given in App. B. 

In order to keep the information on the energy fraction z in each part of the subtraction 
function, the finite part of the real corrections is modified to 



do 



77^4/7 

finite 



1 

2~i 



d$ 



4/7 



1^77-4/7 |2 Q / $ 



4/7 J 



- J2 |Mmb,*j| 2 ©(pi = Zijpi, k = (1 - Zij)pi, {Pk^i}) 



>,j=i 



(4.17) 



It is easily seen that the variable z^, which is defined in Eq. (4.8), plays the role of the 
energy fraction Z{ in the collinear limit for each dipole ij. Again, in the collinear-safe 
case the functions of the subtraction function depend only on the sums Pi + k = pi of 
collinear momenta; in this case we recover Eq. (4.10). 

In the integration of the subtraction function over the photon phase space, we now have 
to leave the integrations over Zij open. The resulting Zij dependence of the integrand is 
most conveniently described with a [...]+ prescription 5 , which separates the soft singularity 
at z^ = 1. The endpoint part at z^ = 1, which results from the full integration over z^, 
exactly corresponds to the contribution of G^ nb \sij) for the collinear-safe case, as given 
in Eq. (4.12), where §ij = IpiPj. The continuum part in involves an integral over 



i(sub)^~ 
Hj 



1 ,1 



5 z ij) 



with 



In m -1 
\ mi ) 



+ {l + z) ln(l -z) + (1-^). 



(4.18) 



The total integrated subtraction part explicitly reads 



/ 



77^4/7 



smg 



+ J\z [g^(~s tJ ,z)] + e( Pl = zp t ,k = (1 - z)p u {p k ^})^. (4.19) 



Owing to the [...]+ prescription, the continuum part is zero if the full integration over z 
is carried out, thereby recovering the collinear-safe case (4.11). 

5 We use the definition dx [f(x)] + g(x) = dx f(x) [g{x) — g(l)]. 
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5 Numerical results 

5.1 Input parameters and setup 

We use the following set of input parameters [34]: 





1.16639 x 10- 5 GeV" 2 , 


a(0) = 


1/137.03599976, 


a s = 


0.1172, 


M w = 


80.423 GeV, 


r w = 


2.118 GeV, 






M z = 


91.1876 GeV, 


r z = 


2.4952 GeV, 






m e = 


0.510998902 x lO" 3 GeV, 


m li = 


0.105658357 GeV, 


m T = 


1.77699 GeV, 


m u = 


0.066 GeV, 


m c = 


1.2 GeV, 


m t = 


174.3 GeV, 


m d = 


0.066 GeV, 


m s = 


0.15 GeV, 


m h = 


4.3 GeV. (5.1) 



If not stated otherwise, the Higgs mass is Mr = 170 GeV. In some cases we alternatively 
use Mh = 130 GeV. The corresponding values for the Higgs-boson decay width Th, which 
have been obtained with the program HDECAY [35], are given by 

T H (M H = 170 GeV) = 0.3834 GeV, T H (M H = 130 GeV) = 0.004995 GeV. (5.2) 

We set the quark- mixing matrix to the unit matrix throughout, but in the limit of massless 
external fermions a non-trivial quark-mixing matrix can be included by a simple rescaling 
of the cross sections. 

Furthermore, we apply a set of recombination and separation cuts: 

1. Bremsstrahlung photons that are closer than 5° to a charged fermion or have less 
energy than 1 GeV are recombined with the charged fermion that is closest in angle. 
This means that in this case before evaluating distributions or applying phase-space 
cuts the momenta of the photon and the fermion are added and associated with the 
fermion, while the photon is discarded. 

2. The following separation cuts are applied to the momenta defined after a possible 
recombination: 

£,>10GeV, #(Z,beam)> 5°, 9(1, 1') > 5°, 9(l,q)>5°, 

E q > 10 GeV, 6(q, beam) > 5°, m(q, q') > 10 GeV, (5.3) 

where an obvious notation for energies E___, angles 9(. . .), and invariant masses 
m(. . .) for leptons / and quarks q is used. 

The separation cuts and input parameters are the same as in Ref. [3] 6 for the processes 
77 — > 4/, so that we reproduce the Born cross sections that we calculated there. In 
particular, we exclude forward and backward scattered charged fermions, because they 
cause collinear singularities. While for final-state quarks these singularities signal a non- 
perturbative regime, for leptons they are in principle cured by finite-mass effects. However, 
we exclude this region by demanding that leptons appear in the detector with finite 

6 There is a misprint in Eq. (6.1) of Rcf. [3]. The value for a s is supposed to be a s = 0.1172 and not 
a s = 1.1172. 
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production angle and energy. Compared to Ref. [10] we use different recombination cuts, 
because, in contrast to e + e~ collisions, the recombination criterion based on invariant 
masses does not lead to collinear-safe observables. This is due to the collinear singularity 
that arises if a charged fermion is collinear to the beam. Even though an appropriate 
cut on the angle between charged fermions and the beam is imposed, it might happen 
that a photon with relatively high energy is recombined with a low-energy fermion that 
is close to the beam. Thus, after recombination, the fermion almost follows the direction 
of the photon and is not affected by the angular cut. Such events are avoided by taking 
a recombination condition based on the angle. 

For the evaluation of the lowest-order matrix elements of 77 — > 4/ and 77 — > 4/7, 
we use the fixed-width scheme, in which the gauge-boson width is introduced in all (i.e. 
time- and space-like) propagators. As argued in Section 2.4 of Ref. [3], this scheme does 
not break gauge invariance for reactions 77 — > 4/ (+7) with massless external fermions. 

The photon spectrum is accounted for by using the parametrization of the program 
CompAZ [4], as described in Section 5 of Ref. [3]. In order to distinguish the cases with 
and without convolution over the photon spectrum, we write ^/s^ and y/s yy for the CM 
energies in these cases, respectively. 

In the numerical integration we generate 2 • 10 7 events for the plots showing the in- 
tegrated cross sections, and 5 • 10 7 events for distributions and for the integrated cross 
sections in Table I. If not stated otherwise, the shown results are based on the subtraction 
method, but have been cross-checked with the slicing approach. Moreover, we have addi- 
tionally checked most of the results by reproducing them within statistical uncertainties 
with our second independent Monte Carlo generator. 

5.2 Integrated cross sections 

In Table 1 we present a survey of integrated cross sections for a leptonic, a hadronic, 
and two semi-leptonic final states, as obtained with the subtraction and slicing methods. 
The cross sections of the semi-leptonic final states differ because of the effective polar- 
izations of the photons resulting from the Compton backscattering (cf. Section 6.3 of 
Ref. [3]). Final states that differ only in the fermion generation (i.e. in their mass values) 
receive the same radiative corrections, since our predictions are based on the massless 
limit for the external fermions and mass singularities cancel after performing a photon 
recombination. The results obtained with the two methods for treating the real correc- 
tions, subtraction ("sub") and slicing ("sli"), are in good agreement. Note that they both 
are implemented in the same Monte Carlo generator, which, thus, yields identical results 
for y/s^j < 170 GeV where the IBA is used. This is the reason why the "sub" and "sli" 
numbers are identical in the case of ^/s^ = 200 GeV with 7 spectrum, where only the 
range y/s^ < 170 GeV is relevant in the convolution. 




In Figure 5 the integrated cross section for 77 — > z/ e e + du including radiative corrections 
is compared with the Born cross section as a function of the CM energy for monochromatic 
photon beams. The "best" curves correspond to the C(a)-corrected cross sections. A 
Higgs boson of Mh = 170 GeV produces a sharp peak in the cross section at */s 77 = 




170 GeV, while for larger energies the corrections are almost independent of the Higgs 
mass. The relative corrections 8 = 0"/o"Born — 1 in the four lower plots of Figure 5 behave 
roughly like the corrections to on-shell W-pair production [15-17]. Close to the W-pair 
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1296.49(52) 


1295.29(62) 


1335.13(38) 


0.09(6) % 



Table 1: Integrated cross sections for different final states and energies with and without 
convolution over the photon spectrum. The third column shows the result obtained with 
the subtraction method and the fourth with the slicing method. The last two columns 
show the Born cross section and the relative difference between subtraction and slicing. 
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Figure 5: Integrated cross section for 77 — > i/ e e + du (the two upper plots) and relative 
radiative corrections (the four lower plots) without convolution over the photon spectrum 
for Higgs masses M H = 130 GeV and 170 GeV. 
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y/s^l GeV 


200 


240 


260 


280 


300 


500 


1000 


TU 


2.0% 


1.9% 


1.3% 


0.8% 


0.7% 


0.5% 


0.5% 



Table 2: Estimates of the TU (5.4) for the C(o;)-corrected cross section of 77 
at various CM energies a/s^T. 



z/ e e + du 



production threshold the corrections are dominated by the Coulomb singularity. For 
higher energies the corrections decrease until they reach about —7% at 1 TeV. In this 
region they are dominated by large logarithms from the Regge and Sudakov domains. 

In Figure 5(c) we also show the comparison with the IBA for a Higgs mass of Mh = 
130 GeV. Since close to the W-pair production threshold the bulk of the corrections is 
due to the Coulomb singularity and since there are no other pronounced corrections, the 
agreement between the two curves is quite good. The very good agreement of the DPA 
and the IBA at ^/s 77 ~ 170 GeV both for semi-leptonic and for hadronic final states (in 
both cases the difference is well below 0.1%) is of course accidental. For the leptonic final 
state the difference is about 0.7%. 

As explained in Section 3.5, the intrinsic uncertainty of the IBA is about ~ ±2%, 
while the DPA accuracy is up to < 0.5% where it is applicable. Since the convolution 
of the hard 77 cross section, in general, involves both the IBA (in the low-energy tail) 
and the DPA (for ^/s 77 > 170 GeV), the uncertainty of our cross-section prediction is in 
the range 0.5—2%, depending on the contribution of the IBA part to the full convolution. 
Denoting the IBA and DPA parts of the full cross section as A<t iba and A<7 DPA (both 
including the corresponding lowest-order contribution, so that A<t IB a + A<t DP a — a), we 
can estimate the theoretical uncertainty (TU) of the corrected cross section a to 

TU = x 2 % + x 0.5%. (5.4) 

a o 

Table 5.2 illustrates this estimate for a few CM energies for 77 — > z/ c e + du. For 
^/s^ < 230 GeV our prediction possesses a TU of ~ 2%, because it is mainly based on the 
IBA, but already for v /s^ > 300 GeV (500 GeV) the IBA contribution is widely suppressed 
so that the DPA uncertainty of < 0.7% (0.5%) sets the precision of our calculation. We 
note, however, that the overall uncertainty of our calculation certainly becomes worse as 
soon as TeV energies for ^/s 77 are dominating because of the relevance of high-energy 
logarithms beyond 0(a). 

In Figure 5(e) the comparison of the full correction with the IBA is shown for a Higgs 
mass of Mh = 170 GeV. The IBA includes the Higgs resonance via an effective coupling 
and reflects the shape of the resonance quite well. 

The cross section including the convolution over the photon spectrum as a function 
of CM energy is shown in Figure 6 for a Higgs mass of M H = 130 GeV and in the lower 
left plot also for Mh = 170 GeV. In the upper plots the integrated cross sections are 
shown, and in the lower plots the corrections relative to the Born cross section. Recall 
that we use the IBA below ^/s 77 = 170 GeV. This means, in particular, that the Higgs 
resonance is calculated from the effective coupling and not from the full DPA in this 
region. The interesting structure in the lower left plot reflects the shape of the photon 
spectrum convoluted with the Higgs resonance. Since the Higgs resonance is very narrow, 
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Figure 6: Integrated cross section for 77 — > z/ c e + du (upper plots) and relative radiative 
corrections (lower plots) including the convolution over the photon spectrum for Higgs 
masses of M H = 130 GeV and 170 GeV (lower left plot). For ^/s^ > 300 GeV (shown on 
the r.h.s.) the "best" curve for Mh = 170 GeV practically coincides with the shown curve 
for M H = 130 GeV. 
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a sizable contribution is only possible if X\X 2 s ee ~ M\ where X\ and x 2 are the energy- 
fractions carried by the photons. The correction is very small at low y^s^T where X\ and 
X2 have to be so large in order to match this condition that the corresponding spectrum 
is extremely small. Increasing ^/s^ allows for lower values of x\ and x 2 . For instance, 
for M H = 130 GeV, the rise at ^/s^ ~ 180 GeV results from a region where both x\ 
and x 2 are in the high-energy tail of the spectrum which is produced by multiple photon 
scattering. The peak at ^/s^ ~ 200 GeV is caused by events where one photon comes from 
the high-energy tail and one from the dominant peak in the photon spectrum. Finally, 
at ^/s^ ^210 GeV both x 1 and x 2 originate from the dominant photon-spectrum peak 
which causes the steep rise until ~ 220 GeV. 



5.3 Distributions 

In Figure 7 we show the invariant-mass distributions for the z/ e e + and du pairs in the 
process 77 — > z/ c e + du, both with and without convolution over the photon spectrum. The 
upper plots show the absolute predictions, and the lower plots the corrections normalized 
to the Born predictions. Since we use y/s^ = 500 GeV or ^/s^ = 500 GeV, the corrections 
are shifted upwards when including the photon spectrum, because the effective energy of 
the photons is lower (cf. Figure 5). The shape of the corrections, however, is hardly 
changed by the convolution over the photon spectrum. As the shape of the corrections 
determine a possible shift of the peak of the invariant-mass distribution, it is of particular 
importance in the determination of the W-boson mass. The measurement of the W-boson 
mass can, e.g., be used for understanding and calibrating the detector of a 77 collider. 

The distribution in the W-boson production angle is sensitive to anomalous couplings. 
In order to set bounds on these couplings it is mandatory to know radiative corrections, 
because both anomalous couplings and radiative corrections typically distort angular dis- 
tributions. The corresponding angular distribution of the du system, which is equal within 
the statistical error to the distribution of the z/ e e + system, is shown in Figure 8. While the 
correction without the photon spectrum is about —9% for W bosons emitted perpendic- 
ular to the beam, the corrections are rather small when including the photon spectrum. 
As already explained above, the cross section is dominated by a region where the 77 CM 
energy is smaller. In fact, the relative correction 5 is accidentally small at ^/s^ ~ 500 GeV 
[cf. Figure 6(d)] and might also become larger if other cuts or event selection procedures 
are applied. 

Figure 9 shows the energy distribution of e + and d for the process 77 — > z/ e e + du. 
The characteristics of the Born cross section, especially the influence of the effective 
polarization of the photons after Compton backscattering, were explained in detail in 
Section 6.3 of Ref. [3]. The relative corrections shown in the lower plots amount to 
a few per cent. For very low and very high energies, where the Born cross section is 
very small, the relative corrections in DPA are not reliable anymore. In this region 
the assumption that doubly-resonant diagrams dominate is not fulfilled. The angular 
distributions for e + and d are shown in Figure 10. The shape of the Born cross section and 
the influence of the photon spectrum were also explained in Section 6.3 of Ref. [3]. Similar 
to the angular distributions of the z/ e e + and du systems, the corrections are maximal in 
a region where the fermions are emitted perpendicular to the beam. However, after 
including the photon spectrum, the corrections almost cancel as can be anticipated from 
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Figure 7: Invariant-mass distributions of the W + and W bosons reconstructed from the 
z/ e e + and du pairs in the process 77 — > i/ e e + du at y/s = 500 GeV. 
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Figure 8: Angular distribution of the W boson reconstructed from the du pair in the 
process 77 — > z/ e e + du at y/s = 500 GeV. 



Figure 6(d) which shows that the corrections to the integrated cross section are almost 
zero at ^fs^ ~ 500 GeV. 

Finally, the energy distribution of the photon in the process 77 — > z/ e e + du+7 is shown 
in Figure 11. The distribution is dominated by the soft-photon pole at k° — > and 
decreases rapidly at higher energies. Comparing the distributions with and without con- 
volution over the photon spectrum, the convolution shifts the curve to lower energies, 
because the initial-state photons already have less energy. 

5.4 Non-collinear-safe observables 

As explained in Section 4.2, the treatment of collinear singularities in non-collinear- 
safe observables deserves some care. Applying the generalizations of the subtraction and 
the slicing methods described above, we now turn to observables without photon recom- 
bination. Apart from that, the same phase-space cuts are applied as before. In Figure 12 
we show the distributions of the ^ e e + , and du pairs in the processes 77 — > z/ e e + du, 

z^/i+du. With photon recombination the leptonic invariant masses of the two processes 
receive the same radiative corrections since the recombination guarantees the necessary 
inclusiveness so that all mass singularities cancel. If the recombination is not applied, 
the distributions change drastically. Note, however, that the recombination is mainly 
a rearrangement of events, and omitting the recombination affects the integrated cross 
section by less than 0.5%. With decreasing invariant masses the relative corrections rise, 
while they are smaller at large invariant masses. The reason is that without recombination 
final-state radiation (which is enhanced by mass logarithms) reduces the invariant mass of 
the reconstructed W boson, thereby shifting events from the dominating resonant region 
to lower invariant mass values. The recombination brings most of these events back to 
the resonance region, because it prevents momentum loss from final-state radiation. The 
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Figure 9: Energy distribution of e + and d in the process 77 — > z/ e e + du at \/i = 500 GeV. 
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Figure 10: Distributions in the production angles of e + and d in the process 77 — > u e e + du 
at V« = 500 GeV. 
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Figure 11: Energy distribution of the final-state photon in the processes 77 — > t / e e + du + 7 
and 77 — > z/ M /i + du + 7 at s/s = 500 GeV. The l.h.s. compares the distributions with and 
without convolution over the photon spectrum (with photon recombination); the r.h.s. 
compares the cases with and without photon recombination (without convolution over 
the photon spectrum). 




Figure 12: Invariant-mass distributions of the W + and W~ boson reconstructed from the 
z/ e e + (z/ At yU + ) pair and du pair in the process 77 — > z/ e e + du (77 — > z/ My u + du) at y/s^ = 
500 GeV, with and without photon recombination. 
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Figure 13: Energy distributions of e + and /i + in the processes 77 — ► z/ e e + du and 77 — > 
z/u// + du at ^/s 77 = 500 GeV, with and without photon recombination. 

l.h.s. of Figure 12 also shows a hierarchy in the mass effects of the outgoing leptons as 
the slope for the v e e + pair is much steeper than the slope for the z/^/i + pair due to the 
smaller mass of e + . The plot on the r.h.s. shows that the corrections for the du pair are 
not as large as for the pair on the l.h.s., because the remaining mass terms behave 
like Qjlnmj, where Qf denotes the charge of the fermion /. We also note that the cor- 
rections are smallest in the case with photon recombination because of the cancellation 
of all mass singularities. 

The photon recombination also affects the energy distributions of the fermions. Fig- 
ure 13 shows this distribution for e + and fi + in the processes 77 — > z/ e e + du and 77 — > 
z/^/i+du with and without recombination. In the former case the curves coincide, as ex- 
plained above. The recombination attributes the photon to a fermion and, thus, shifts 
events to higher energies of the fermion. The mass-singular effect, which appears without 
recombination, is again larger for e + than for 

The effect of the photon recombination on the photon-energy spectrum is shown in 
Figure 11. Without recombination the distribution is shifted to higher photon energies 
because the recombination transfers events to the bin with zero photon energy. The 
difference is again bigger for the process 77 — > i/ e e + du than for 77 — > z/^/i + du, since the 
mass-singular logarithms of e + are larger. 

6 Summary 

In this paper we have described a calculation of the 0{a) electroweak radiative cor- 
rections to 77 — > WW — > 4/ in the electroweak Standard Model in the double-pole ap- 
proximation (DPA). Technically, we follow the strategy of the RacoonWW Monte Carlo 
event generator for the corresponding e + e~ reaction. This means, virtual corrections are 
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treated in DPA and are decomposed into factorizable and non-factorizable contributions, 
while real-photonic corrections are based on complete lowest-order matrix elements for 
77 — > 4/7. The combination of virtual and real corrections is done in two different ways: 
by using the dipole subtraction method or by applying phase-space slicing. 

A detailed survey of numerical results for the 0(a) corrections has been given, com- 
prising results on integrated cross sections as well as angular, energy, and W-invariant- 
mass distributions. In the W-pair threshold region the corrections are dominated by the 
Coulomb singularity and are, thus, positive and of the order of a few per cent. For in- 
creasing 77 scattering energies the corrections become more and more negative and reach 
about —10% in the TeV range for integrated cross sections. For large scattering angles, 
where the Born cross section is relatively small, the impact of the corrections is usually 
larger. Since the convolution with realistic photon beam spectra effectively reduces the 
hard scattering energy, the size of the corrections is usually somewhat reduced compared 
to the situation with monochromatic photon beams. Typically, collinear-safe observables 
(i.e. where mass-singular logarithms cancel due to sufficient inclusiveness) receive correc- 
tions of a few per cent for energies of the e~e~ system before Compton backscattering up 
to 1 TeV. As expected, non-collinear-safe observables receive very large corrections (tens 
of per cent) because of the existence of logarithmic mass singularities. 

The radiative corrections are implemented in a Monte Carlo generator called 
COFFER77, which optionally includes anomalous triple and quartic gauge-boson cou- 
plings in addition and performs a convolution over realistic spectra of the photon beams. 
The construction of this generator and lowest-order results obtained with it have already 
been described in a previous publication [3]. At present, C0FFER77 is the only event gen- 
erator that includes both the decays of the W bosons and radiative corrections, thereby 
defining the state-of-the-art in the description of the processes 77 — > WW — > 4/(+7). 



Appendix 

A Transformation of the coefficient functions Fj 

In this appendix we describe the transformation of the coefficient functions Fj for the 
factorizable virtual corrections (3.9) that transforms all Fj into the helicity amplitudes of 
the on-shell process 77 — * WW. 

The 36 SME M] yW+W ~ of Ref. [16], which fix the coefficient functions Fj by Eq. (3.6), 
are defined for 36 different helicity configurations which can be enumerated with a single 
index /, 

My w+w ~(X 1 ,\2;X + ,X-) = M jl , Z = (A 1 ,A 2 ,A+,A_), (A.l) 

where j, I — 1, . . . , 36. The 36x36 matrix M is explicitly obtained by inserting momenta 
and polarization vectors into the 36 independent SME _A/fJ 7W+w of the 83 structures 
defined in Eqs. (5)-(9) of Ref. [16]. 
If we transform the Fj according to 

36 

F l = Y l F j M jh (A.2) 
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the function F[ is the helicity amplitude for the on-shell process 77 —> WW corresponding 
to the helicity configuration I = (X 1 , X 2 , A + , A_). As such, it can be well approximated by 
the generalized Fourier series described in Section 3.2.2. It is important to notice that in 
Ref. [16] the scattering plane spanned by the beam axes and the produced W bosons was 
rotated into the (x 1 , x 3 )-plane, so that the SME J\4^ w+W depend only on s and cos 6, 
or equivalently on s and t. Since, thus, the matrix M is a function of s and t, also the 
new functions Fi depend only on s and t, but not on the azimuthal angle of the scattering 
plane or other on kinematical variables. According to Eq. (A. 2), the SME M.j transform 
as 

36 

Mi = ^(M- 1 )^, (A.3) 
i=i 

where M~ l denotes the inverse matrix of M. By construction, the transformation de- 
couples the different helicity channels of 77 — > WW. When including the W decays in 
the SME, as done in Eq. (A.3), this decoupling is somewhat disguised for the W-boson 
polarizations, but still valid for the photon helicities. This means that the new SME 
Mi consist of four subsets, each of which contributes only for one of the four different 
polarization combinations (Ai,A 2 ) of the photons. In practice, we have evaluated and 
simplified the matrix M and the new SME Mi analytically as much as possible. 

B Dipole subtraction for non-collinear-safe photonic final-state radiation 

In Section 4.2.2 we have collected the relevant formulas for the generalization of the 
dipole subtraction method to non-collinear-safe observables. Here we describe the details 
of their derivation. Specifically, we focus on the situation of light charged particles in the 
final state only; the more general case of massive particles and of charged particles in the 
initial state will be worked out elsewhere [36]. Although not made explicit in the main 
text, we keep track of the polarizations of the outgoing particles. 

(i) Subtraction of singularities 

Generically the schematic form of the subtraction procedure to integrate the squared 
matrix element X)a 7 I-^i| 2 (summed over the photon polarizations A 7 ) for real photon 
radiation over the (N + l)-particle phase space d<3>i reads 

J d$! ]T \M,\ 2 = J d$x \M,\ 2 - \M suh A + jd<§o ® (/ [dk] \M suh \ 2 ^) , (A.l) 

where d$o is a phase-space element of the corresponding non-radiative process and [dk] 
includes the photonic phase space that leads to the soft and collinear singularities. The 
sign "®" indicates that this factorization, in general, is not an ordinary product, but 
may contain also summations and convolutions. The two contributions involving the 
subtraction function |.M su b| 2 have to cancel each other, however, they will be evaluated 
separately. The subtraction function is constructed in such a way that the difference 
(Xa 7 \-Mi\ 2 — |-M S ub| 2 ) can be safely integrated over d<3>i numerically and that the sin- 
gular integration of |-M S ub| 2 over [dk] can be carried out analytically, followed by a safe 
numerical integration over d^o- 
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In the dipole subtraction formalism, the subtraction function is given by [22,23] 



|AW$ i)| 2 = -£ QiOiQjOjJg^pipuPj, k) \M (%,f, | 2 , (A.2) 

where the sum runs over all emitter-spectator pairs ij, which are called dipoles. Recall 
that both i and j are final-state particles in our case. The relative charges are denoted Q / 
(/ = hj)i an d the sign factors 07 = ±1 correspond to the charge flow (07 = +1 for anti- 
fermions, 07 = —1 for fermions). The summation over r = ±1 accounts for a possible flip 
in the helicity Kj of the emitter %. The singular behaviour of the subtraction function is 
contained in the radiator functions gf^r\piiPji k), which depend on the emitter, spectator, 
and photon momenta Pi, Pj, and k, respectively. In the limit of small fermion masses 
the functions gfj^f 1 are related to the function gf^ of Eq. (4.7) for the unpolarized 
case by gfjf^ = g^°\ gfj^ — 0- The squared lowest-order matrix element |.Mo| 2 of 
the corresponding non-radiative process enters the subtraction function with modified 
emitter and spectator momenta pi, pj, as defined in Eq. (4.9). The momenta are related 
by pi+pj + k = pi+pj, where all the other particle momenta pk and pk, entering |.Mi| 2 and 
|.Mo| 2 , respectively, are the same, pk = Pk- The modified momenta are constructed in such 
a way that pi — > pi + k in the collinear limit (pik — > 0). Since we deal with light external 
fermions only, we set all masses mj of external fermions to zero whenever possible. This 
means that m/ = can be consistently used in the integral / d$i (j2\~ 
but that the readded contribution J[dk] |-M su b| 2 in general contains mass-singular terms 
of the form alum/. 

In collinear-safe observables, and only those are considered for light fermions in 
Refs. [22,23], a collinear fermion-photon system is treated as one quasi-particle, i.e., in 
the limit where a charged fermion i and 7 become collinear only the sum Pi + k enters the 
procedures of implementing phase-space selection cuts or of sorting an event into a his- 
togram bin of a differential distribution. Technically this level of inclusiveness is reached 
by photon recombination, a procedure that assigns the photon to the nearest charged 
particle if it is close enough to it. Of course, different variants for such an algorithm 
are possible, similar to jet algorithms in QCD. The recombination guarantees that for 
each photon radiation cone around a charged particle % the energy fraction of Eq. (4.14), 
z i = Pi/iPi + 1S fully integrated over. According to the KLN theorem [33], no mass 
singularity connected with final-state radiation remains. Collinear safety facilitates the 
actual application of the subtraction procedure as indicated in Eq. (A.l). In this case 
the events resulting from the contributions of |A4 su b| 2 can be consistently regarded as 
A^-particle final states of the non-radiative process with particle momenta as going into 
M.q (^cMj) , i-e. the emitter and spectator momenta are given by pi, pj, respectively. 

Owing to pi — > Pi + k in the collinear limits, the difference (Da 7 I-^i| 2 — |A^ su b| 2 ) can be 
integrated over all collinear regions, because all events that differ only in the value of 
enter cuts or histograms in the same way. The implicit full integration over all Zi in the 
collinear cones, on the other hand, implies that in the analytical integration of |.M su b| 2 
over [dk] the Zi integrations can be carried out over the whole Zi range. 

In non-collinear-safe observables, not all photons within arbitrarily narrow collinear 
cones around outgoing charged particles are treated inclusively. For a fixed cone axis 
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the integration over the corresponding variable is constrained by a phase-space cut or 
by the boundary of a histogram bin. Consequently, mass- singular contributions of the 
form alnrrii remain in the integral. Technically this means that the information on the 
variables Z{ has to be exploited in the subtraction procedure of Eq. (A.l). The variable 
that takes over the role of Zi in the individual dipole contributions in |-M su b| 2 is as 
defined in Eq. (4.8), because in the collinear limit it behaves Thus, the integral 

/ d$i (Jlx., \M-i\ 2 — |.Msub| 2 ) can be performed over the whole phase space if the events 
associated with |.M su b| 2 are treated as (N + l)-particle events with momenta Pi = Zijpi, 
Pj = pj, and k — (1 — Zij)pi. This modification, in turn, requires a generalization in 
the evaluation of the second subtraction term on the r.h.s. of Eq. (A.l), because now the 
integral over z; L j implicitly contained in [dk] depends on the cuts that define the observable. 

(ii) Integration of singular contributions 

For a final-state emitter % and a final-state spectator j with masses TOj and rrij the 
integral of (PiiPj, k) over [dk] is proportional to 

4> } (4) = ^1= r^-y^ l ZM ^ii9tr\Pi^k), (A.3) 
2JXij J yi Jzi{yij) 

where Pfj = (pi + pj) 2 and the definitions of Section 4.1 of Ref. [22] are used. There the 

results for G < f^\P 2 ^) with generic or light masses are given in Eqs. (4.10) and (3.7), re- 
spectively. In order to leave the integration over zij open, the order of the two integrations 
has to be interchanged, and the integral solely taken over yij is needed, 



*7t(l%- %) = -3= dy« (1 - y i3 ) g^\ Pi , Pj , k). (A.4) 

9. / A._ Jyi(zij) 



Note that the function Qf^\pfj, z) is not needed for finite photon mass A, because the 
soft singularity appearing at z — > 1 can be split off by employing a [...]+ prescription in 
the variable z, 



S§?\P&,z) = G^\P 2 )5(l - z) + [Gt™\P^z)\ + . (A.5) 

This procedure shifts the soft singularity into the quantity Gf^\Pf 3 ), which is already 
known from Ref. [22]. Moreover, the generalization to non-collinear-safe integrals simply 
reduces to the extra term Qf^\Pf^z) , which cancels out for collinear-safe integrals 
where the full z- integration is carried out. 

In the limit m ; — > and for rrij = A = the boundary of the integration is given 

by 

y^) = mK i' z \ »W = i. (A.6) 

and the functions relevant for the integrand gf^ behave as (see Section 4.1 of Ref. [22]) 

Pf. 

Pik = -y2/«, R ij(y) = 1 ~y, r ij(y) = i- ( A -7) 



?(sub) / 73,2 
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The evaluation of Eq. (A. 4) becomes very simple and yields 



In 



2 

mf 



- 1 



+(l + z) ln(l-z), 



(A.8) 



where Pff{y) is the splitting function of Eq. (4.16). Equation (A.8) is correct up to terms 
suppressed by factors of Wj. For completeness, we repeat the form of the full integral 
Gf^\Pfy in the case of light masses, 



7T 



Tf(sub)/ 



1 

2' 



with the auxiliary function C of Eq. (4.13). The results for the functions Qf"^ and G\ 
for the unpolarized case, as given in Eqs. (4.18) and (4.12), are obtained by summing over 
the variable r = ±1 which accounts for the spin-flip, 



(A.9) 



-((sub) 



^j(sub) _ (sub) , Msnh) 



/-y(sub) ^ysuo) x-y 



y(sub) 



(sub) 



(A.10) 



Finally, we give the explicit form of the ij contribution \M. su b,ij(®i)\ 2 to the phase- 
space integral of the subtraction function, 



a 



-^QiViQjVj I d&o,ij\Mo(pi,pj;TKi^- 



d$i |M>ub,ij($i)| 

x jcgf (J3)e(So«) + jT dz [Qt^\p^ z)] + e( Pl = z p u k = (i - z)p t , {p k ^}) 

(A.n) 

The arguments of the step functions 0(. . .) indicate on which momenta phase-space cuts 
are imposed. 
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